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! Abstract 



The unknown parameters of the recently-proposed {Int J. Geom. Meth. Mod. Phys. 1, 751 

[2004]) Castro-Mahecha model of the imaginary parts (Aj) of the nontrivial Riemann zeros are 

r~| . the phases (a^) and the frequency parameter (7) of the Weierstrass function oi fractal dimension 

-4— > ■ 

^ ' D = I and the turning points {xj) of the supersymmetric potential-squared $^(a;) — which in- 

corporates the smooth Wu-Sprung potential {Ptiys. Rev. E 48, 2595 [1993]), giving the average 

■ level density of the Riemann zeros. We conduct numerical investigations to estimate/determine 
00 ! 

00 ■ these parameters — as well as a parameter (cj) we introduce to scale the fractal contribution. Our 

primary analyses involve two sets of coupled equations: one set being of the form $^(xj) = \j, and 

■ the other set corresponding to the fractal extension — according to an ansatz of Castro and Ma- 

(—j ■ hecha - of the Comtet-Bandrauk-Campbell (CBC) quasi-classical quantization conditions for good 

' 2 1/2 

g I supersymmetry, p(^^/2) f-x ^{x -x'^-^^/^ ~ analyses suggest the possibility strongly 

^ ' that 7 converges to its theoretical lower bound of 1, and the possibility that all the phases {ak) 

^ ' should be set to zero. We also uncover interesting formulas for certain fractal turning points. 

^; 
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I. INTRODUCTION 



The precise nature of the zeros of the Riemann zeta function is of paramount mathe- 
matical [1], as well as physical interest P]- Of course, perhaps the most famous, not yet 
fully resolved, mathematical conjecture (the Riemann Hypothesis P|) is that all the non- 
trivial zeros {'r]j + iXj) — the trivial zeros simply equalling —2j, j = 1,2,3, .. . — have real 
components rji = 

Castro and Mahecha recently proposed a supersymmetric implementation of the Wu- 
Sprung model of the A/s, employing a Weierstrass (continuous and nowhere differentiable) 



function "6] for the fractal structure of dimension -D = f that Wu and Sprung had found for 
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their local one-dimensional potential (V) (cf. 0,0)- The Wu-Sprung potential V — which 
generates the smooth average level density obeyed by the Riemann zeros — satisfied Abel's 
integral equation 0, eq. (6)], and was written imphcitly as 0, eq. (7)], 

Here Vq = 3.10073^ ^ 9.74123. 

The Weierstrass fractal function employed by Castro and Mahecha J, eq. (81)], 



1 _ 1X7*= 

W{x, 7, D,a,) = o^^^^e^™'' , (2) 

has both an unknown (countably) infinite set of phases 27rak and an unknown parameter 
7, which determines the frequencies 7^. Further, D is the fractal dimension, which Castro 
and Mahecha — following the box-counting argument of Wu and Sprung [5] (cf. |9j) — 
took to be |. (Castro and Mahecha also noted ^, sec. 6] that this choice of D = | does, 
appealingly, yield [the omnipresent] j noise.) The variable x parameterizes the inverted 
Wu-Sprung potential Vws{x), which in accordance with the supersymmetric requirement 
that the energy be zero in the ground state, was translated to be zero at a; = 0. 

Rather than directly addressing the formidable problem "of factoring the ordinary 
Schrodinger equation studied by Wu-Sprung", Castro and Mahecha (CM) postulated the 



model 







|,eq. (82)], 

<^\x) = Vwsix) + ^[W{x, D, 7, ak) + W{-x, D, 7, a^) + c.c] + 0o, (3) 

with 

$2(x,) = A,. (4) 



(c.c. denotes complex conjugation.) The additive constant 0o = —Vo is chosen so that, in 
accordance with supersymmetric principles, $^(0) = 0. The "turning points" xj (where the 
potential energy is equal to the total energy), along with the phases and the parameter 
7, are the unknown (countably infinite) parameters of the CM model, which we hope here 
to estimate/determine, by satisfying the equations 0, simultaneously with the quantization 
conditions ^ eq. (87)], 

^.(^., A,) = 2-- 1^ ^-J^-^^iyq^dx' = jn. (5) 

Castro and Mahecha proposed (cf. [lO|, eq. (36)]) that these relations would constitute the 
fractal extension of the fermionic phase path integral approximation (the CBC [Comtet- 
Bandrauk-Campbell] formula (cf. [o[ eq. (1)] B eg. (5)])) — the supersymmetric counter- 
part of the well-known WKB quantization rule UJ, eq. (2)]. Here = = f and j runs 
over the positive integers. 

The Castro-Mahecha model Q was cast in a fractal supersymmetric quantum-mechanical 
(SUSY-QM) setting, in an effort to implement the Hilbert-Polya proposal to prove the 
Riemann Hypothesis, under which an Hermitian operator, hypothetically, can be found to 
reproduce the X/s as its spectrum. It had appeared to CM quite difficult to directly solve 
the SUSY Schrodinger equation [J, eq. (84)], 

+ $)(_i?(/3) + $)^(+)(a;) = A,^]+)(x), (6) 
where h = 2m = 1, Aj-^"* = A^- the ifjj^^ are eigenfunctions, and the fractional derivative 

V^''^F(t) = - -T -^^^dt' (7) 

i_. Tii-p)dtJ_^it-tr^ 

is intended (cf. |lO|). (M. Trott suggested that one might directly pursue the solution of 
this system, most effectively, using a Griinwald-Letnikov discrete approximation ^^.) CM, 
therefore, had recourse to the CBC formula (jSj) as an approximation to the implementation 
of the quantum inverse scattering method, which would yield the potential $. Given an 
exact potential, one could construct the SUSY Schrodinger equation, giving the Aj's for its 
eigenvalues. This would realize — CM claimed — the Hilbert-Polya proposal to prove the 
Riemann Hypothesis. 

CM reasoned that if the potential energy has the shape of a fractal curve, then it is very 
natural to assume that the kinetic energy operator should be a fractal; namely it should be 
expressed in terms of fractal derivatives. 



In sec. m we begin our numerical analyses by studying the solutions to Q- Then, in 
sec. mil we extend the analyses to include (simultaneously) the two sets of equations Q 
and (0). The fractal Wu-Sprung potential V is examined in sec. I1V[ in particular its turning 
points, for a number of which interesting formulas are obtained. In sec. |3 we expand the 
CM model to include a certain scaling factor. A discussion of our analyses is contained in 
sec. IVII We also include an Appendix outlining recent work of D. Dominici concerning the 
inversion of the Wu-Sprung potential Q| . 



II. PRELIMINARY ANALYSES 

A. Turning points xj^^ — for smooth potential — held fixed 

Our initial numerical analyses of the CM model consisted of, for various integral cut-offs 
m on the running index k in the Weierstrass fractal function summation Q, finding an 
optimizing set of a^s, while keeping the x/s fixed and equal to {n of) the turning points for 
the translated inverted smooth Wu-Sprung potential, that is, those x^-^'*'s for which (cf. 0), 

Vws{xf) + (po = Xj- (8) 

The set of estimated a^s and 7 would be chosen, so as to minimize the sum of squared 
deviations — since it did not appear possible to achieve complete equality — between the 
LHS (left-hand side) and RHS of (jH)). Computationally-speaking we employed the NMinimize 
command of Mathematica, using the "DifferentialEvolution" and "InitialPoints" options, 
and for the initial points only increasing sets of values for the a^'s, in order to reduce search 
times. Also, the necessary inversion of the Wu-Sprung potential Q — in order to obtain 
Vws{x) — was performed using the Interpolation command of Mathematica (cf. Appendix). 

We, first, found certain evidence, at least in this initial context, that suggested taking 
7 = 3 = 2D. (In 4], CM wrote that "it would be intriguing to see if 7 = 1 -|- ^ 1.618, the 
inverse of the golden mean, since the golden mean appears in the theory of quantum noise 
related to the Riemann Hypothesis" (cf. ^, p. 481]). 

Initially, we tried the case n = 100, m = 15, insisting that the phases be ordered from 
lowest to highest. Our results were encouraging, and we proceeded to successively higher 
integral values of m, reaching m = 100. (We make the obvious comment that here the 
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FIG. 1: The one hundred ordered phases (divided by 2tt) of the Weierstrass fractal function 
obtained by minimizing the sum-of-squares deviation from the imaginary parts of the first one 
hundred nontrivial Riemann zeros of the Castro-Mahecha supersymmetric model 

number m of unknowns [a^] equals the number of [nonlinear] equations.) For that case, 
the best fit to the first 100 Aj's yielded a remarkably small sum-of-squares deviations of 
2.68927 ■ 10~^^. The first seventy-seven of the hundred ordered phases (that is, 2nak, k = 
1, . . . , 70) were all equal to ^ = 27r/3 to, typically, seven decimal places (Fig. H}. When we 
inserted the turning points (xj^^) derived from the (smooth) inverted Wu-Sprung potential 
for the bounds of integration — that is, those satisfying (jHJ — into the fractal CBC formula 
(0), we found that these quantization relations were clearly not fully met. (In general, the 
CBC quantization conditions appear to hold for a number of "shape-invariant" potentials 
[ul (cf. Q).) In Fig. 121 we plot the (CBC) ratio of the computed value of the real part 
of the LHS of ^ to its RHS (that is, jir). We see that for the lower values of j, the LHS's 
exceed jir, while for the higher values, they are less than jn — all in an apparently smooth 
monotonically decreasing manner. 

We, then, were able to extend our computer analyses to the case n = m = 200, with 
very analogous results (Fig. (The sum of squared deviations was quite small, as in the 
n = m = 100 analysis, 3.48792 ■ 10~^^.) Most of the estimated values of at were once again 
I to high precision, this time for = 1, . . . , 144. In Fig. E] we present the analogue of Fig. |21 
for the n = m = 200 case. 

In regard to the phases of the Weierstrass fractal function, Castro and Mahecha sec. 
5] had written that it "is warranted to see if the statistical distribution of these phases ak has 
any bearing to random matrix theory (the circular unitary random matrix ensemble) and 
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FIG. 2: Ratios of the (real parts of the) computed values of the left-hand side of the fractal CBC 
quantization conditions © — based on the (fixed) turning points, those Xj for which Xj = Vwsixj), 
of the smooth inverted (translated) Wu-Sprung potential and the scaled phases shown in Fig. ^ — 
to the right-hand side of © (that is, jvr). The ratios would all be unity if the conditions were fully 
satisfied. 
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FIG. 3: The two hundred ordered phases (divided by 2tt) of the Weierstrass fractal function (jJJ 
obtained by minimizing the sum-of-squares deviation from the imaginary parts of the first two 
hundred nontrivial Riemann zeros of the Castro-Mahecha supersymmetric model ^ 

the recent studies of quantum phase- locking, entanglement, Ramanujan sums and cyclotomy 
studied by H". 

When we proceeded to the still higher-order case n = 300, m = 300, we found superficially, 
at least, different types of results. In six different runs (using different random seeds) of our 
Mathematica program, we obtained sum-of-squares deviations no greater than 
1.79962-10^^^, all six yielding essentially the same-looking plot. The smallest sum-of-squares 
was 1.62071 ■ 10"^'^, corresponding to Fig. |S| The first seventy-four scaled phases (a^) were 
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FIG. 4: Ratios, for the case n = 200, of the (real parts of the) computed values of the left-hand 
side of the fractal CBC quantization conditions © to the right-hand side of © (that is, jvr). The 
ratios would all be unity if the conditions were fully satisfied. 
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FIG. 5: The three hundred ordered phases (divided by 27r) of the Weierstrass fractal function 
obtained by minimizing the sum-of-squares deviation from the imaginary parts of the first three 
hundred nontrivial Riemann zeros of the Castro-Mahecha supersymmetric model 

all equal to \ (rather than |, as previously) to high accuracy. Let us note, in any case 
though, that for either ay, = | or |, the factors e^'^*"'' in the expansion (j2I) of the Weierstrass 
fractal function have zero real part. 



B. Comment 



What it appeared, in retrospect, we had actually principally accomplished above (Figs.Hl 
El) was to find sets of phases (27rafc) for various n (100, 200 and 300), such that the real 
part of the Weierstrass fractal function Q is very close to zero — as revealed by plots — 
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for any x (or 7 > 1), be it a turning point or not. So, these results immediately pertain to 
the Weierstrass fractal function and not the Riemann zeros themselves. Nevertheless, our 
estimated phases (a/c's) can be employed as the first step of an iterative procedure — as 
discussed in the next section. 

Since the Castro-Mahecha supersymmetric model of the Riemann zeros requires the si- 
multaneous obtaining/fitting of phases (a^'s), the parameter 7, and turning points (x/s), 
we turn to this problem now. 



III. JOINT ESTIMATION OF PHASES (a/,) AND TURNING POINTS {xj) 
A. Iterative approach 



1. n = m = 100 analyses 

We would like now to regard the two sets of equations ^ — the one set we have tried to 
satisfy in the above analyses — and (jSj) as a coupled system and attempt to simultaneously 
solve for the phases {ak) and turning points (xj) (still holding 7 = 3). In this context, we 
again considered the case n = 100, m = 100, and took the solution portrayed in Fig. [T]as 
the first step of an iterative procedure. Using this set of a^.^'^'s, we searched for a new set of 
turning points {x^^\j = 1, . . . , 100), such that the CBC quantization conditions (0) would 
all he exactly (up to our numerical accuracy) satisfied. Then, with this revised set of Xj's, 
we hoped to find anew the set of a^^^^s satisfying Q. By reiterating this two-step procedure 
a sufficiently large number of times, we anticipated possibly arriving at a final set of a^f^^'s 
and a;^-^'''s simultaneously satisfying (jH} and ^ (cf. A test for convergence to such a 

final solution would be that the two sets of parameters would only change negligibly upon 
further iterations. 

At the very beginning of the iterative procedure we have just outlined, we sought to 

(2) 

obtain a (hypothetical, second) set of one hundred turning points {xj ), which would render 
the CBC quantization conditions ^ completely satisfied when the first set of phases (27ra^^'') 
— those depicted (in their scaled form) in Fig. Q — were employed. However, this proved to 
be not totally doable. For example, for the original value of x^qq = 15.315, the corresponding 

CBC ratio was 0.926293. But when we tried to adjust x^qq to increase this ratio closer to the 

(2) 

ideal value of 1, we were only able to obtain a s%/it improvement by going to x^qo ~ 14.3452, 



CBC ratio 




FIG. 6: CBC ratio — incorporating the phases shown in Fig. ^ — as a function of the hundreth 
turning point (xioo)- The original value of x^^qq, giving a ratio of 0.926293, was 15.315. This can 

(2) 

be (only) slightly improved to 0.937452, by choosing x^qq = 14.3452. 
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FIG. 7: The original 100 turning points {xj ), forming the higher /dominant line, and the 100 

(2) 

adjusted/revised turning points {xj ) that yield CBC ratios optimally closer to unity — given the 
phases shown in Fig. ^ 

for which the ratio was 0.937452. We show this phenomenon in Fig. 

In all one hundred cases, when we sought to drive the CBC ratios closer to unity, the 
new turning points chosen, x^^\ were smaller than the originals x^-^"* (Fig. Ej). (There is some 

(2) (2) 

anomalous behavior apparent, in that x\-{ < Xjq .) For the first seventy-one cases, the new 
turning points {xj ) do yield the desired CBC ratio of 1 (Fig.lHj). (We have not been able to 
obtain a revised set of a^^^'s based on the 100 a;^-^'''s to closely satisfy Q, as we had with the 
100 xj^'^'s.) When the revisedset of turning points xj^"* was employed, along with the original 
set of phases a^l^^s (Fig. [TJ, instead of highly accurate predictions of the A/s, we obtained 
signficant discrepancies. (The sum-of-squares measure of fit to (jHI) jumped enormously from 
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FIG. 8: CBC ratio as a function of the first turning point (x^^"*). The value of x^^\ which gave 



(Fig. a ratio of 3.48049, was 1.30083. This can be optimally improved to 1, by choosing x\ 
0.141784. 
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FIG. 9: The difference, divided by jvr, of the j-th Riemann zero minus that value for the zero 
predicted from © using the revised set of turning points {xj ) (Fig. [TJ and the original phases 

4' 
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2.68927- 10"^"^ to 129157.) In Fig. 121 we show the differences (all positive) of the X/s minus 
their predicted values, divided by jn. 

In this paper, we have not adhered to the apparent misstatement in jj, eq. (89)], that 
the Riemann zeros (Aj) should be "properly normalized" ("unfolded" [3, eq. (6)]) as 

A,-^logA„ (9) 

but rather conducted our analyses in terms of the "raw" unnormalized Riemann zeros them- 
selves (A/s), similarly to the work of Wu and Sprung y (cf. Q). 
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B. Non-iterative (simultaneous) approach 

1. n = m = 8 analyses 

We commenced an additional series of analyses — in lieu of the iterative scheme just 
proposed and discussed (but not fully pursued, due to apparent computational demands). 
For the case n = m = 8, we tried to simultaneously estimate the scaled phases (a^) and 
turning points [xj), as well as the frequency parameter 7 (for a total of seventeen param- 
eters). We sought to minimize — as we will throughout the remainder of this study — a 
certain sum-of-squares, one of the two addends being composed of the squared differences 
between the LHS and RHS of the (SUSY-QM) equations (jH)) and the other addend of the 
squared differences between the LHS and RHS of the (CBC) equations (0). (We simply, in 
effect, weighted the two sums equally, but there certainly appears — in retrospect — to be 
opportunity for refinement in the choice of weights, since the second (CBC) sum seemed 
to consistently contribute considerably more weight.) We also — as mentioned — allowed 
the parameter 7 that determines the frequencies to vary (between 1 and 5), rather than 
holding it fixed at 3 (as we had preliminarily thought there was some numerical evidence 
for doing so earlier (sec. ITT|)). 

The minimum sum-of-squares we achieved was 181.6414. A predominant part of this sum 
— 130.236 — was attributable to the lack of fit to the CBC quantization conditions ©. 
Remarkably, all eight phases (a^) were estimated to be essentially zero (or, effectively, the 
same 27r). The largest deviation of these eight estimates from or 2t[ was, in fact, only 
8.347 ■ 10^^. (We had abandoned also here our earlier computational constraint that the 
phases be monotonically nondecreasing and only required, numerically, that they lie between 
and 27r.) For 7 we obtained an estimate of 1.4119. The eight estimated turning points {xj) 
were (0.662734, 1.17022, 1.58516, 2.28689, 2.32938, 3.18343, 3.21295, 3.30561), which can be 
compared with (1.30083, 1.87866, 2.20626, 2.64243, 2.84142, 3.20489, 3.4613, 3.64459), based 
on the smooth Wu-Sprung potential (jH)). Based on these two sets of turning points — and 
ak = 0, k = 1, . . . ,8 and 7 = 1.41119 — the corresponding CBC ratios (in the same order of 
presentation) are (1.80006, 1.55891, 1.3248, 1.32142, 1.22064, 1.23564, 1.23128, 1.21166) and 
(1.46964, 1.34903, 1.19796, 1.26407, 1.14572, 1.23298, 1.201, 1.16905). We can see that the 
latter set of CBC ratios is clearly superior/more desirable — giving rise to a sum-of-squares 
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deviation of only 83.5937 vs. the 130.236 mentioned-before — but use of the original set 
of turning points (xj^'') yields a much larger sum-of-squares for deviations from the other 
(SUSY-QM) set of equations © of 620.878, as opposed to the 51.4054 = 181.6414 - 130.236 
we obtained. 

2. n = m = 7 analyses 

We had also conducted some parallel n = m = 7 analyses. We obtained from Mathemat- 
ica, three minima (with figures of merit of 69.6285, 140.887 and 171.939). The corresponding 
estimates of 7 were 1.16585, 1.34436 and 1.4099. Of the 21 (7 x 3) phases estimated, none 
was greater than 0.000210 removed from either or 27r. 

We will find, apparently importantly, in Sec. IV Al that for this n = m = 7 case, we can 
obtain a much improved figure of merit of 13.703 by incorporating one additional param- 
eter (cr) into the SUSY-QM model (jHI), one which scales the Weierstrass fractal function 
iy(x,7, Before doing so, however, we will examine, in detail, the fractal Wu-Sprung 

potential. 

IV. ANALYSES OF THE WU-SPRUNG FRACTAL POTENTIAL 

We implemented, for the case n = 100, with an algorithm provided by M. Trott, the 
dressing transformation for the inversion of the eigenvalues (Aj's), employed by van Zyl and 
Hutchinson jsljin order to obtain the fractal form VwSfrad^) Wu-Sprung potential 

(Fig. HOD (cf. U Fig. 2], Figs. 1, 3]). In Fig. ITD we show the non-smooth component 
(the "residual" in statistical parlance) of the fractal WS potential. In other words, we have 
subtracted away from FigHUl the smooth WS potential Vwsi^)- We have tried to manipulate 
the Weierstrass fractal function as incorporated into ^ — assuming all ak = — so 
as to closely resemble Fig. ^2 A simply visual, non-rigorous analysis lead us to present 
Fig. El in which we take 7 = 2.3 and multiply the Weierstrass fractal function term in (jH)) 
by 5 and subtract 10. (Possibly, the CM model might be modified, somewhat in line with 
this observation. So, it would appear, that even though the Weierstrass fractal function, 
with its fractal dimension D set to | and the fractal component of the VwSfrac share the 
same fractal dimension, they seem to differ considerably in certain of their overall [scale. 
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location] properties. We will explore this issue in Sec. IVl) 



Vwsf„c (x) 




FIG. 10: Fractal Wu-Sprung potential, VwSfraS^)^ based on the first 100 Riemann zeros 
A. Turning points {xj^"''^) of the fractal Wu-Sprung potential 

Using the results of the dressing transformation, we initially considered taking for the 
values of Xj^"''^, those closest to the original smooth Wu-Sprung turning point (x^^), for 
which Xj = VwSfraci^j^'^'^)^ new/revised turning point to employ in the CM model. 

Pertaining to the issue, for the edification of the reader, we will note that in an e-mail 
exchange of ours with M. Trott, he had commented: "But they (van Zyl and Hutchinson 
use an exact solution of the Schrodinger equation. So in their approach they never 
need turning points. This bring us back to my remark from the beginning of this project. 
What is a 'turning point' in a fractal potential? If you take the smallest XMin, such that 



Vws£„c (x) -Vws (x) 




FIG. 11: Non-smooth component of the Wu-Sprung fractal potential VwSfraS-^)^ that is, 
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5 Wfrac(x)-10 
25r 




FIG. 12: Weierstrass fractal function contribution to ^ multiplied by 5, then diminished by 10, 
with 7 = 2.3 and = 0, with intention of adjusting the Weierstrass fractal function to resemble 
Fig. ^2 as closely as possible 

Jl^" - V{x) = nh with e = V{xMin), then, because of the fractal nature, the barrier 
could be very small and the tunnelling contributions from larger x could be quite significant." 
In response, Castro first wrote:" Since the Weierstrass function oscillates, there are many 
possible choices for the a^J^"^, so you have to choose for the turning points those which are 
the closest to the ones [x^j \ obtained from the smooth Wu and Sprung potential alone. 
This process will determine the x^^'^'^ points uniquely". Trott further suggested: "From 
an ordinary quantum mechanics point of view using the 'nearest' ones seems unconvincing 
to me. Already for the simple case of two wells separated by a wall, the classical JWKB 
formulas get nontrivial corrections" (cf. (However, in retrospect, it now appears that 

the most effective manner in which to resolve this issue — at least, in pragmatic terms — 
is to choose that x^^'^'^ for which the associated fractal CBC relation (0) is most closely 
satisfied.) 

We, in fact, found nine possible values of ^I"'' ^^at is (0.0359371, 0.224349, 0.403363, 
0.60336, 0.75138, 0.902281, 0.949646, 1.5688, 1.63925), satisfying the turning point relation 
VV5;,,,(4r') = -^1- The correponding CBC ratios (Fig. ^ — (0., 0.766538, 1.04323, 
0.871811, 1.42656, 1.26132, 1.30892, 1.09161, 1.27263) seemed, except obviously for the 
first, reasonably well-behaved. For the second Riemann zero (A2), the ten turning points 
found were (0.432856, 0.57922, 0.787248, 0.848629, 1.13689, 1.19066, 1.54042, 1.66097, 
1.88178, 2.32703) with corresponding CBC ratios (Fig. ED (0.966061, 0.848087, 1.2481, 
1.18487, 1.24222, 1.26174, 1.14487, 1.35361, 1.24891, 1.16204). For A3, the sixteen turning 
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FIG. 13: CBC ratios based on the fractal Wu-Sprung potential VwSfrad^)- "^^^ i^ii^e {h = 1, . . . ,9) 
turning points (xf*^""^) for which VwSfraci^i^'^'^) ~ indicated. 
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FIG. 14: CBC ratios based on the fractal Wu-Sprung potential VwSfraS^)- '^^^ ^^^^ = !)•••! 10) 
turning points (x2^"'^) for which VwSfraS^il'^'^) ~ "^^ indicated. 

points (0.450109, 0.563705, 0.816881, 1.0183, 1.10108, 1.22306, 1.38131, 1.39341, 1.5228, 
1.67159, 1.86394, 1.90039, 1.90039, 2.30223, 2.35184, 2.79924) with corresponding CBC ra- 
tios (Fig.Cni) of (0.765244, 0.689671, 1.01581, 1.15142, 1.09523, 1.15534, 1.08533, 1.08337, 
1.04846, 1.22877, 1.14428, 1.15193, 1.15193, 1.07091, 1.0915, 1.02537). We also show results 
(Fig.Cni) for x^or'' ^i^h turning points (9.11597, 9.22462, 9.4841, 9.9511, 10.019) and fractal 
CBC ratios (0.980756, 0.971956, 0.97427, 0.948974, 0.94763). 

B. Formulas for certain fractal turning points 

Continuing along these lines, we found that particular value (0.949646) of x for 
which VwSfraci^) ~ '^i) which, of the nine solutions, was, at the same time. 
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FIG. 15: CBC ratios based on the fractal Wu-Sprung potential VwSjraX'^)- "^^^ sixteen {h = 
1, . . . , 16) turning points (xg^""^) for which Vw s fraS^i]^'^) ~ indicated. 
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FIG. 16: CBC ratios based on the fractal Wu-Sprung potential VwSfraS^)- '^^^ (/i = 1, . . . , 5) 
turning points {x{q'^) for which VwSfraci^io^) ~ -^3 indicated. 

smaller than x^^^ = 1.30083, and closest to it. With the use of "Plouffe's inverter" 
( http: / / pi.lacim.uqam.ca/eng/) , we obtained a remarkable formula that would appear to 
apply to it: 

xf™'^ = 0.949646 ^ 0.99999996 ■ 10"^^, (10) 
where t ~ 1.8392867 is the tribonacci constant and GS is the G elf ond- Schneider constant, 

Proceeding further along these lines, the value {x = 1.660974) we obtained for x^^"''^, that 
satisfied VwSfrad^) — -^2, and was smallest and closest to x''2^ = 1.87866 was expressible as 

g2A3 

x^J''"' = 1.660974 ^ 1.0000028005 ■ 10"^^ ^ ^. (11) 

' Uog'(2 + v^) ^ ^ 
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Additionally, 



STrott 

= 1.9003895 ^ 0.99999998 ■ 10"^°e^^^ (12) 

' (2 + ^3)2' ^ ^ 



where Trott = 0.010841015122311136 is the Trott constant, and = 2.20626. (The Trott 
constant has the property that it is invariant, in a certain sense. If one expands the number 



digit by digit as if it were a continued fraction, then the number remains the same 
Further, 

e4A7 

= 2.3843247 ^ 1.00000003037 ■ lO"^^.^ — (13) 

t^Cahen^\og({5) 

where Cahen ~ 0.6434105462883 is the Cahen constant, used in the theory of continued 
fractions. (0:4^^ = 2.64243 and ( is the Riemann zeta function — exact formulas for which 
are known for even integral arguments.) 
Also, 

(A/o\ 4v^ 5A9 

= 2.8338417 ^ 1.00000004769 • 10'^°^ ^^ ' . (14) 

■'^(12) 

{xl^""^ is rather close to x^^^ = 2.84142. ({3) is known to be irrational and is sometimes 
referred to as "Apery's constant" 122^.) 



V. SCALING THE WEIERSTRASS FRACTAL FUNCTION 

The analyses underlying Figs. ^2 and ^1 and their comparison, indicated to us that pos- 
sibly the CM model of the Riemann zeros might be enhanced by multiplying the Weierstrass 
fractal contribution by some (new) scaling factor (a). (Since in Sec. m we were essen- 
tially finding parameters to set the Weierstrass function to zero, the issue of scaling was not 
germane there.) So, we modified the type of analyses reported in Sees. 1111 B II and 1111 B 2l to 
estimate a, as well as the (scaled) phases (a^), the turning points (xj) and the frequency 
parameter (7). 



A. n = m = 7 analyses 



We use exactly the same criterion of best (unsealed) fit, employed previously, that is the 
total sum-of-squares between the LHS's and RHS's of ^ and (0). The sum-of-squares we 
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obtained was 13.703 (8.48114 coming from the lack of fit to the CBC quantization conditions 
(0)), much superior to the best fit of 69.6285 in Sec. HUB 21 This, of course, lends strong 
support for the relevance of incorporating the scaling parameter a, which was set to 3.92036, 
along with 7 = 2.18081 — both estimates being quite consistent with Fig. ^1 The seven 
estimated phases were (0.915274, 6.28319, 6.28319, 1.20429, 5.33637, 0.700917, 0.0) — three 
of them being essentially or 27r 6.28318531. The seven estimated turning points were 
(0.321253, 0.676572, 0.936234, 1.65921, 1.79613, 2.1688, 2.18378), considerably smaller than 
in Sees. HUB 21 and HUB II The associated CBC ratios — all now relatively close to unity 
— are (1.1386, 1.16732, 0.951104, 1.13392, 1.02881, 1.067, 1.06951). 

At this stage, we could have continued this n = m = 7 analysis to see if we could obtain 
further minima with figures of merit less than 13.703. But now knowing that Mathematica 
could produce results for = m = 7, we were interested in seeing how far we could "push" 
Mathematica for cases n = m > 7. So, rather than now dwelling on a particular n = 
m case, we kept seeking higher-order instances. (At this point, we treated each analysis 
independently, that is we did not employ results of earlier [lower-dimensional] analyses as 
initial guesses.) We did obtain single minima for n = m = 8 and n = m = 9, but the figures 
of merit were rather weak (39.1008 and 85.689), so we do not detail the results here. 

B. n = m = W analyses 

Now we expand to the case of ten phases and ten turning points. The (relatively 
quite excellent) sum-of-squares we obtained was 2.86609 (2.68872 coming from the lack 
of fit to the CBC quantization conditions (0)). The estimate of 7 was 1.30466 and of a, 
1.49575. The ten estimated phases were (6.28319, 6.28319, 1.27142, 0.0493542, 6.28319, 
0.0555025, 0.0000178313, 6.28319, 0., 6.00153 • lO^^S), most of them being essentially or 
2tt ^ 6.28318531. The ten estimated turning points were (0.38402, 0.665384, 0.915631, 
1.32033, 1.42245, 1.83676, 1.95767, 2.33568, 2.68742, 3.06041), considerably smaller than 
in Sees. IIII B 21 and IIII B 11 where the scaling parameter (a) was absent. The associated 
CBC ratios — most now relatively close to unity — are (1.27306, 1.142, 1.00359, 1.03838, 
0.969627, 1.02007, 1.0058, 0.978438, 1.00323, 0.984601). 

So, an important unresolved question is whether or not ideally all the phases (a^) should 
be set to zero, or whether the fact that all the phases were estimated as zero in Sec. HUB II 
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was principally a manifestation of the need for a scaling parameter greater than unity. This 
is because if one sets all the phases to zero, the term-by-term interference in the Weierstrass 
fractal function is minimized and its overall contribution/infiuence/magnitude increased. (I 
thank Carlos Castro for this last observation.) 

1. Zero phases 

To address this last-discussed issue, still within the n — m — 10 scenario, we a priori set 
all ten phases (ak) to zero, leaving us with twelve parameters (ten turning points, plus 7 
and a) to estimate. We obtained four minima, one relatively inferior (24.0743) and three, all 
yielding approximately 7.2358 (6.758 of this from the CBC equations) — but all more than 
the 2.86609 obtained immediately above — for the sum-of-squares, and all with essentially 
the same set of estimates. The estimates were 7 = 1.25274, a = 1.23289 (both interest- 
ingly close to one another) and for the turning points (0.406011, 0.72083, 0.88735, 1.21503, 
1.38215, 1.72066, 2.07865, 2.24067, 2.81678, 2.91862). The CBC ratios were (1.4716, 1.20207, 
1.04624, 1.0322, 0.953197, 0.985181, 0.987027, 0.9602, 1.03323, 0.995425). 

So, no yet fully convincing evidence for the proposition that all the phases (ak) should 
be set to zero, has appeared in this n = m = 10 suite of analyses (though Castro believes 
this is an interesting hypothesis well worth pursuing). 

C. n = m = 12 

We were able to obtain two minima {independently-generated with different random seeds) 
for this case, having associated sums-of-squares of 4.7436 (4.51953 attributable to the fit to 
the CBC quantization conditions) and 4.88101 (4.59739 so-attributable). 

These two minima had rather similiar estimates of 7, a, as well as of the twelve turning 
points. (For the smaller of the two, 7 and a were estimated as 1.16902 and 1.90025, and for 
the larer, 1.17209 and 1.79825.) The turning points for the smaller (of the two small) minima 
were: (0.370638, 0.518695, 0.65557, 1.3734, 1.49501, 1.90065, 2.13024, 2.25476, 2.83058, 
2.96561, 3.18434, 3.69962), and for the larger of the two (0.371207, 0.526829, 0.697696, 
1.33586, 1.45628, 1.91474, 2.13455, 2.25897, 2.76173, 2.88164, 3.10636, 3.50429). The phases 
(27rafc) for the smaller of the two were: (5.06388, 1.19176, 5.10283, 5.10239, 1.54197, 4.87505, 
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5.75582, 5.07329, 6.28319, 6.28319, 0.491806, 0.) (three being essentially or 27r) and 
(1.36199, 6.28319, 5.32029, 3.14014, 6.28319, 0.932408, 1.48183, 6.2407, 0.811536, 6.28186, 
0.0160344, 0.) (five now being essentially or quite close to or 27r). The CBC ratios based on 
the smallest (of the three n = m = 10) minima were: (1.45746, 1.09132, 0.891351, 1.03399, 
0.967626, 1.0088, 1.01311, 0.981312, 1.01104, 0.986514, 1.00474, 1.00054). 

D. n = m = 15 

Mathematica yielded a minimum having a sum-of-squares of 10.8781, 10.1905 of the total 
stemming from deviations from the idealized CBC relationships. 

The estimate of 7 was 1.10999 and of the scaling parameter cr, 2.39239. The fifteen 
phases were estimated to be (5.85115, 1.90972, 6.07707, 3.38707, 5.79917, 3.01828, 5.21112, 
0.972974, 6.28319, 4.06528, 5.13804, 1.09178, 5.42993, 0.0000400221, 0.887475) and the fif- 
teen turning points, (0.427238, 0.575873, 0.674962, 0.964599, 1.67025, 1.86543, 1.99707, 
2.09904, 2.87102, 2.96694, 3.12156, 3.26476, 3.37853, 3.42232, 3.66246). The resultant 
CBC ratios were (1.62295, 1.205, 0.970261, 0.919054, 0.964421, 1.03902, 1.01992, 0.974149, 
1.00741, 0.984101, 1.00373, 1.01898, 1.01317, 0.977452, 0.996592). 

Subsequently, we obtained a second minimum (starting with different random seeds), with 
a slightly inferior sum-of-squares of 11.6345 (10.9979 due to failure to completely satisfy the 
CBC relations). The estimate (1.10638) of 7 obtained was notably similar to the first esti- 
mate (1.10999), while the estimate of a was 1.92895 vs. 2.39239, previously. The estimated 
phases were (1.44074, 6.28318, 1.66813, 0.542437, 2.10471, 4.76708, 2.12488, 0.0381873, 
1.27475, 0.672527, 1.10953, 1.92716, 5.96083, 0.00013761, 0.) and the estimated turning 
points, (0.447782, 0.605638, 0.711853, 1.03604, 1.65431, 1.87966, 2.02289, 2.13206, 2.85323, 
2.95601, 3.11569, 3.25825, 3.37216, 3.41606, 3.70782). The correlation coefficient between 
the two sets of turning points was extremely high, 0.999796, while that between the two sets 
of phases was positive, but apparently quite weak, 0.174294. (The sets of turning points are 
constrained to be strictly increasing, while the phases are not.) This can be (maximally) 
improved to 0.471313 by adding 6 = n to the second set of phases fFig. ll7|) . The correlation 
between the two sets of CBC ratios was extremely high, 0.997186. (Obviously, no ordering 
constraints has been placed upon them.) 
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FIG. 17: Correlation coefficients between tlie first set of pfiases for the n = m = 15 case, and the 
second set of phases, with the addition to its members of an angle G [0, 27r]. The highest plateau 
of level 0.471313 extends from = 1.56108995 to 4.1583035 

E. n = m = 20 

Initially, we obtained a minimum with a somewhat disappointingly large (in comparison 
with previous analyses) sum-of-squares of 60.1355 (57.0834 stemming from the failure to 
fully satisfy the CBC relations). The estimate of 7 was 1.03784. Then, we obtained a 
somewhat superior sum-of-squares of 53.6492 (50.1484 coming from the lack of complete fit 
to the CBC relations). The estimate of 7, 1.049, was quite similar. 

Ultimately though, we found a much superior solution with a sum-of-squares equal to 
11.4018, of which 10.3952 was attributable to the failure to fully fit the CBC relations. The 
estimate of 7 was 1.08285 and of a, 1.23124. The twenty phases were (6.28319, 5.12918, 
1.82172, 0.627236, 2.28354, 6.28057, 2.09377, 2.27574, 6.28319, 6.28319, 1.39899, 1.05003, 
4.08259, 0.00871371, 6.28319, 5.20743, 0.00527643, 6.28319, 3.08389 • 10"^ 0.00224054) — 
half of them being quite proximate to or 27r — and the twenty turning points, (0.430261, 
0.581042, 0.682297, 0.990577, 1.66733, 1.8869, 2.03975, 2.16318, 2.7836, 2.8999, 3.11796, 
3.36483, 3.61917, 3.71294, 4.1214, 4.23216, 4.39019, 4.59062, 5.39774, 5.49934). (Possibly, 
the relatively large number of phases near either or 27r is a computational artifact, since 
these are the extreme values between which the phases are constrained to lie.) Addition- 
ally, the resultant CBC ratios were (1.63002, 1.21105, 0.97568, 0.925948, 0.959536, 1.03755, 
1.02111, 0.978159, 1.01279, 0.983753, 0.998834, 1.01199, 1.00774, 0.97948, 1.01524, 1.00143, 
0.997587, 0.993286, 1.00448, 0.994628). 

We have been investigating the n = m = 25 case (both with unconstrained phases and 
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phases set to zero) — with some programming improvements provided by M. Trott — but 
do not have any specific resuhs to report presently. (It appears that several days on a 
PowerMac are required to locate a single relative minimum.) 



1. Decreasing trend in estimates of j 

We see an undeniable trend in the several analyses above. As the number of unknowns 
{n + m + 2) increases, the parameter 7 decreases, possibly indicative of a convergence to its 
theoretical lower bound for the Weierstrass fractal function of 1. (For n = m = 7, our one 
estimate of 7 was 2.18081; for n = m = 10, the estimate was 1.30466; for n = m = 12, we 
had two estimates of 1.16902 and 1.17209; forn = m = 15, 1.10999 and 1.10638; and for 
n = m = 20, 1.03784, 1.049 and 1.08285.) Szulga and Molz have shown that in the limit 
7 — > 1 (with random phases, in the case D = ^ before us here), "the Mandelbrot- Weierstrass 
p.oce. . a complex .ac.o„a, B.own,a„ .o^on" Q (cf. Q). In t.e Ca... 

Mahecha model, one takes the real part of the Mandelbrot- Weierstrass function (which CM 
term the Weierstrass function). To statistically test whether or not the estimated series of 
phases we obtain is random in nature, one can apply "Rao's spacing test" |28]. (Also, the 
Rayleigh test and the Kuiper's V test are sometimes employed.) Applying this test to the 
first two sets of twenty phases reported for our n = m = 20 analyses in Sec. IV El we obtained 
statistics of 192.199 and 158.86, respectively. Any result greater than 192.17 is statistically 
significant at the p = 0.001 level of probability, and any result greater than 154.31, at the 
p = 0.10 level js^j Table II]. There is, thus, rather strong evidence that the hypothesis that 
the phases are random should be rejected. 

VI. DISCUSSION 

On the scientifically important question of possible falsifiability, the Castro-Mahecha 



supersymmetric model of the nontrivial Riemann zeros — and that of Wu and Sprung 
— would fail (as Carlos Castro indicated) if it could be demonstrated that the shape of the 
potential were multifractal, so that the fractal dimension would no longer be a constant {D = 
|). We are not cognizant of any evidence of this, though. However, it might be of some value 
to formally test for such a possibility. (Though we are aware of no multifractal counterpart 
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FIG. 18: Absolute value (|$^q(x)|) of the wavefunction corresponding to the tenth eigenvalue Aio 
and the tenth turning point xio = 2.96694 reported for the n = m = 15 analysis (smaller of the 
two minima) in Sec. IV PI 



to the Weierstrass function, Castro suggested that one could sum (or possibly integrate) 
over the fractal dimension (D).) Additionally, still within a "unifractal" framework, Castro 
put forth the idea of integrating (thus, generalizing) the Weierstrass fractal function — with 
fractal dimension D held constant at | — over the parameter 7, using some appropriate 
weighting function f{j). 

The presumptive Hilbert-Polya operator would have eigenvalues equal to the Aj's. It 
might be of interest to attempt to estimate its eigenfunctions within the same postulated 
quasi-classical (CBC) framework employed in the CM model, and in our analyses above. 
In Sec. 6.3 of [2^ (also eq. (56)]) there are presented formulas for the quasi-classical 
wave functions. The appropriate fractal extensions of these formulas — following the ansatz 
employed in (jSI) — would need to be employed (cf. jsD])- An initial effort of ours along 
these lines — by simple way of a single illustration — yielded Fig. ^| In the context of 

zeros and a fractal potential", 

Wu and Sprung noted that "All the wave functions decay exponentially beyond the turning 
point, and thus quickly die away. This implies that the lower energy levels have very small 
influence on the potential beyong their turning points . . . Clearly, the potential in the low 
energy range has greater 'responsibility' than that in a higher energy range, thus it has more 
structure . . . Based on this argument, the finest structure will be determined by the wave 
number of the last wave function." 

It would be of interest to estimate all the n wavefunctions in our various n = m analyses 
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FIG. 19: Counterpart to Fig. based simply on the smooth Wu-Sprung potential, without any 
fractal contribution, and original (non-fractal) formulas for the quasi-classical eigenfunctions 

above, and examine whether or not the fractal SUSY Schrodinger equation (jHI) is (at least, 
approximately) satisfied. This would be, in some sense, a test of the appropriateness of the 
fractal extension of the CBC relations postulated by Castro and Mahecha (adopted by them 
in order to avoid having to directly solve ©). (It would also, of course, be of interest to 
attempt to solve © more directly, rather than having to resort to the fractal CBC ansatz 
of CM.) 

In Fig. El we show the counterpart to Fig. ^1 based simply on the smooth Wu-Sprung 
potential,_without any fractal contribution, and original (non- fractal) formulas (|23, eq. 
(6.40)], jnl, eq. (56)]) for the quasi-classical eigenfunctions. It would be an interesting 
exercise to further duplicate some of the analyses conducted here, but with the replacement 
of the fractal CBC equations (0) with the (conventional) CBC equations. 

Throughout the minimization analyses above, for the sake of simplicity (and to avoid the 
possible confusion of issues), we have relied upon a quite elementary least-squares fitting. 
Modifying this approach to give more equable weighting to the two sets of coupled equations, 
and perhaps to lower and higher Riemann zeros and turning points might prove useful. 

CM believed "that their fractal SUSY QM model, once the optimum value for the ampli- 
tude factor 7 is known, has a great chance of truly reproducing the zeta zeros, and proving 
the [Riemann Hypothesis], by simply establishing a one-to-one correspondence among the 
values of the infinite phases of our Weierstrass function with the zeta zeros" jj, sec. 6]. We 
have tried to make a contribution in such a direction here. 

Castro suggested that the most salient /interesting aspect of our several analyses above 
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had emerged in Sec. IV B 11 where for the case n = m = 10 we had set the ten phases to zero 
(doing so which had been strongly suggested by the resuhs of the immediately preceding 
analysis in which the phases were not constrained, and allowed to vary between and 
2n, but found, nevertheless, to be predominantly or 27r). The estimates of the frequency 
parameter 7 and the scaling parameter cr were quite intriguingly close (1.25274 and 1.23289). 
However, the figure of merit, 7.2358, in this zero-phase analysis was somewhat inferior to 
that [2.86609] in that unconstrained n = m = 10 analysis, which rather induced us here not 
to subsequently pursue solely further analyses {n = m > 10) in which the phases had been 
a priori set to zero — but rather let them freely vary. Also, we suspected, as previously 
mentioned, that the predominance of O's or 27r's might be, in some way, a computational 
artifact, since these are the imposed limits on the phases. (The minimization procedure, we 
speculated, might drive the phases to these end points, beyond which it could not further 
proceed to seek improvements, and would therefore terminate there. Castro suggested that 
setting the range of possible phases to be [— vr, tt] would obviate this particular problem — 
but possibly introduce new ones. We intend to investigate such a direction.) 

Perhaps somewhat relatedly, as one of their "spectral speculations". Berry and Keating 



32L p. 260] remarked that the "Maslov phases associated with the [classical periodic orbits 



of the Riemann dynamics] are also peculiar: they are all tt. The result appears paradoxical 
in view of the relations between these phases and the winding numbers of the stable and 
unstable manifolds associated with periodic orbits , but finds an explanation in a scheme of 
Connes" . (The Maslov phases appear when one modifies the Bohr-Sommerfield quantization 
rules by equating the WKB integrals (orbits) to jn plus a Maslov phase factor, instead of 

JTT.) 

Castro suggested that the possibility the Weierstrass fractal function might assume the 
form of classical Brownian motion was an unappealing one, seeing that this would apparently 
imply that the phases (a^) must be random. (In Sec. IV Ell we deduced certain statistical 
evidence that the estimated phases [allowed to freely vary] were not, in fact, random. Of 
course, we know this to be strictly true, since the phases are constructed by our estimation 
procedure.) 

Appendix 

At our suggestion, D. Dominici investigated the problem of inverting the Wu-Sprung po- 
tential [5, eq. (7)]. His analyses have been presented in a "preliminary report" (l4| . 
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FIG. 20: Inverted smooth Wu-Sprung potential — using the Mathematica Interpolation com- 
mand on a data table of resolution ^ — and the approximation (dashed lines) to the inverted 
potential using the first three terms (|17|) of the Dominici series expansion (|15|) 



Dominici gives explicitly the first ten coefficients (ofc) of a power series expansion, 

Vwsix) = Vo + J:Zia,{7Tx)V-\-Voy-\ (15) 



where 



(16) 



(For the critical value Vb = 27r, above which Vws{x) is single-valued, the expansion (fT^ is 
no longer valid. A separate analysis [l^ Sec. 2.2] is, then, of interest.) Using just the first 
three of these coefficients, 



tti = U] 02 



8 9 28 o 

Is" + -9-" • 



(17) 



we obtain Fig. 1201 

In the analyses reported in the body of this paper, we had relied upon the Interpolation 
command of Mathematica for the purpose of inversion, but we anticipate employing the 
results of Dominici in our further work. C. Krattenthaler has indicated that for the 
purpose of studying the asymptotic behavior of the coefficients of the power series, it would 
be best to employ the methodology of Lagrange inversion. 
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